Stress-strain relations in bulk metallic glasses and colloidal dispersions 
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A comparison is made between the nonlinear rheological response of bulk metallic glass formers 
and of colloidal dispersions. Stress-strain curves measured after switch-on of constant deformation 
rates are analyzed quantitatively using a schematic model of mode coupling theory generalized to 
homogeneous and incompressible flows. A mapping between metallic and dispersion rheology is 
possible when stresses are rescaled by an entropic scale, accumulated strains by geometrical factors, 
and rates by the intrinsic relaxation time. Exploiting this similarity and the possibility to directly 
observe individual colloidal particles, we investigate the structural distortions in the colloidal system 
using confocal microscopy. The distortions exhibit the (from elasticity theory) expected quadrupolar 
but also a strong isotropic component. 

PACS numbers: 64.70.P-, 61.25.Mv, 82.70.Dd, 47.50.Cd 



The nonlinear rheological response of melts and com- 
plex fluids directly relates to the internal motion of the 
constituents. Flow causes alEne motion on average, but 
also affects the fluctuations of particle positions and mo- 
menta. This non-affine motion changes the diffusive 
properties of the particles [2[, affects the stress fluc- 
tuations and distributions [3|, and generally determines 
the macroscopic response of the material [J]. The lat- 
ter is strongly affected because solid structures, either 
stable or metastable are broken up and the system is flu- 
idized. Rheology is thus perfectly suited to investigate 
glass forming materials. 

Bulk metallic glasses exhibit unique mechanical prop- 
erties of technological importance, and have therefore 
been the focus of intense investigations [H . Their micro- 
scopic deformation modes, however, remain incompletely 
understood While metallic melts deform plastically 
on the atomic scale, which determines their casting prop- 
erties, in the arrested state, they typically fail as brittle 
materials by shear banding. 

The intrinsic entities of deformation in amorphous 
metals under flow are still unknown. Their investigation 
is, however, possible using colloidal dispersions. This 
opens the way for experiments such as (confocal) mi- 
croscopy [l|, 0, 0, @1 , echo-diffusing wave spectroscopy @ , 
and scattering techniques [l3l to observe the particle mo- 
tions and arrangements under flow. The present letter 
first establishes a close link between metal and disper- 
sion rheology by a quantitative combined modeling of 
their transient stress response. The rheological response 
is then related to the strain-induced structural response 
by observing the local shear-induced distortions of the 
colloid packing via microscopy and simulations. 



Theoretical progress on glass rheology often proceeds 
via considering specific choices of local particle moves 
and deducing the macroscopic response from them. Free 
volume and shear-transformation-zone theory are two 
prominent approaches [ll|, |12| , and recent progress con- 
cerns long-ranged couplings between individual excita- 



tions |13j . Generally, excitations of quadrupolar angular 
dependence are considered based on arguments from con- 
tinuum elasticity theory. Mode coupling theory (MOT), 
on the other hand, addresses the problem starting from a 
microscopic description of the quiescent structural relax- 
ation in a supercooled liquid |14l | , generalizing it to flow- 
driven systems [l5| . The description relies on (low-order) 
correlation functions which capture the average behav- 
ior of unspecified deformation entities. Non-Markovian 
effects play an important role, although they are only 
added occasionally in other approaches, e.g. to capture 
spectral broadening [l^ . The spatial structures involved 
in the wavevector dependent memory kernels of MCT 
reveal a localization of particles in cooperatively formed 
cages. The localization length quantifies Lindemann's in- 
sight that in a solid just at its limit of stability particles 
rattle around their average sites by about a tenth of the 
average spacing fv^ . 

The present study is aimed at metallic and colloidal flu- 
ids close to their glass transitions, where non-Markovian 
contributions to the dynamics are most pronounced. We 
investigate whether the averaged correlation-function de- 
scription of the effects of plastic deformations on struc- 
tural relaxation can rationalize the behaviors of two 
widely different kinds of liquids. While the dispersion 
is close to its MCT glass transition, the metal alloy is 
below the of MCT [ll . Our work thus addresses the 
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question about the generic flow effects in an averaged cor- 
relation function description of the structural relaxation. 

We start from the generalized Green-Kubo relation for 
the transient (deviatoric) stress tensor cr(t) [l^. For an 
incompressible, homogeneous, and constant strain rate 
tensor it becomes UM- 



The schematic model is completed by an equation of 
motion for the density correlator [/]), which con- 
tains a memory kernel describing non-Markovian friction. 
It models the feedback mechanism in supercooled liq- 
uids that structural rearrangements slow down because 
of higher friction, which itself grows because of slower 
structural rearrangements. This leads to the idealized 



cr(t, 7) = / dt' dB{t - t') /dt git -~t',[B]) . (1) glass transition of MCT [IJ . We use the schematic 



p(7) 



It yields the resulting stress after switch-on of flow. The 
Finger tensor B describes the affinc straining of particle 
separations. The stress correlation function g{t, [B]) cap- 
tures the response of the material. Equation ([T]) entails 
the principles of locality, causality, and material objectiv- 
ity (the rheological response is independent of arbitrary, 
including time-dependent, global rotations) [Til. [20|. 

MCT-ITT captures the internal and non-affine parti- 
cle motion using transient density correlation functions. 
Because of the strong local coupling of the structural dy- 
namics [1^, a single correlation function $(^, [I]) can be 
used. The resulting schematic model [21| states that 
transient stress fluctuations decay by particle interac- 
tions which leads to a quadratic relation 



git, [B])^v,iIit))'S>'it, [I]) 



(2) 



Strain enters the model via the invariants of the Fin- 
ger tensor, which combine to the accumulated strain /(t) 
[l9j. Incomprcssibility leads to a quadratic expansion 
lit) ~ ijt)"^ + . . . with a nominal strain rate 7, which 
suffices in the cases considered here, where the steady 
states are reached for strains around 0.3. (Two differ- 
ent deformation geometries will be considered: simple 
shear and uniaxial compression, resulting in 7 = 7s and 
7 = y/Sju, respectively.) The stress-density coupling ver- 
tex Vail it)) gives the elasticity stored in the structure. 
Affine deformations destroy it gradually, which suggests 
that Vailit)) decreases with increasing strain f2]| : 



Va (/(t)) = < (1 - Iit)/il) cxp [-/(i)/7* 



(3) 



The stress-scale is set by w*. In 'repulsive glass form- 
ers', where the harsh repulsive interactions dominate the 
231 . microscopic MCT finds it to be of en- 
12 l l. For hard sphere systems at the glass 
viz. the thermal en- 
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tropic nature 

transition, v* scales with ksT/R^, 
ergy per particle volume. Equation ([3]) includes two char- 
acteristic strain values. They have been calculated for 
hard disks in two dimensions [2lj, but serve as fitting 
parameters in the schematic model. The value of 7* in- 
dicates when stress fluctuations become anti-correlated, 
gi[I] > 7*) < m. The strain 7** corresponds to the 
limit over which stresses can remain correlated in a glass. 
In fluid states, plastic decay in $(t) cuts off the transient 
stress relaxation earlier. The two strains 7*, 7** are ma- 
terial constants connected to the melting/ cage-breaking 
process (l4. 17|. 



model |2l| , which includes the suppression of memory by 
flow and has been used frequently. 




FIG. 1. Stress vs strain curves of Vitreloy 1 (cr = ctxx, sym- 
bols) for different rates, from Ref. [U, with schematic MCT 
fits (lines; parameters in Tab. the legend in (a) gives the 
compression rates 7„ in s~^. Panel (a) shows data at the glass 
transition T ^ Tg = 62SK, and panel (b) at T = 613K. 

In the following, we first show that the model describes 
compressional stresses a^x measured on the bulk metallic 
glass former Zr4i.2Tii3.8Cui2.5NiioBe22.5 (Vitreloy 1) in 
the vicinity of its glass transition temperature Tg [2J]. 
Then we use the model to map the Vitreloy data to shear 
stress-strain curves measured in colloidal dispersions of 
thermosensitive PNiPAM microgels immersed in water, 
at packing fraction cj) ~ 0.60 close to their colloidal glass 
transition |25| . Finally, we investigate the strain- induced 



changes in the local structure by observing the motion of 
the colloids by confocal microscopy. 

Figure [T] shows stress-strain relations of Vitreloy 1 for 



two temperatures at and below Tg « 623K [2J]. At first. 



the stress increases linearly with accumulated strain, 
axx = 2Goo7u^, which allows the easy identification of 
the elastic shear modulus Goo- According to MCT, stress 
and Goo are scaled in units of ksT / R?, with R an ap- 
propriate particle radius. We chose the covalent radius 
of the majority element, i.e. i?2r=17'5pm. This gives 
Goo ~ ^^ksT / R^ , which describes an increase of the 
stress up to order unity in rescaled units. The flow rate 
becomes relevant when the Weissenberg number Pe= 7r 
indicates that the inverse rate and the internal relax- 
ation time become comparable. We estimate the latter 
from the model fits by the Maxwell relation r = 770 /Goo 
with the Newtonian viscosity 770. For Pe> 1, a clear 
maximum in stress develops (Fig. [5]) because stress fiuc- 
tuations become anti-correlated for larger strains. The 
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stress-correlation function g{t, [B]) in Eq. ([2]) becomes 
negative, and the integral, which gives the stress accord- 
ing to Eq. ([TJ, decreases. The maximum lies at strain 7^,. 
Increasing Pe, the overshoot increases and in some curves 
a shallow minimum develops before at large strains the 
stress approaches its value on the flow-curve. The min- 
imum is absent in the schematic model, where 7** gives 
the strain when the flow-curve asymptote is monotoni- 



cally reached [21 
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r = 613K 


162 


0.0894 


0.122 


T = 623K 


160 


0.0762 


0.115 


Colloids 




(7.) 




(ji = 0.60 


90 


0.10 


0.129 



TABLE I. Parameters in Eq. Q for the fits to the metal- 
lic glass Vitreloy 1 and colloidal dispersion PNiPAM; mean 
values given of the strains; see [2^ for Fjj^ parameters. 

According to MCT, the rescaled stresses in metallic 
melts and colloidal solutions should obey similar rela- 
tions with accumulated strain even though the absolute 
values differ by about nine orders of magnitude [13, HI] . 
Rescaled hy ksT / B? , with i?2r= 175pm in Vitreloy 1 and 
i?H=90nm for the colloidal microgels, data from the two 
systems can be compared quantitatively. This is shown in 
Fig. [21 where the rescaled stresses are shown as a function 
of time rescaled by the strain rate 7. The different ge- 
ometry, uniaxial for the metallic versus shear for the col- 
loidal fluid, enters the relation between the true and the 
nominal strain rate : 7 = 7s and 7 = •\/37m- The char- 
acteristic strain of the position of the stress-overshoot 
for both systems then almost agrees. This implies, fol- 
lowing Lindemann's insight, that the local cages which 
freeze-in at the structural glass transition are compara- 
ble in soft and hard systems, if rescaled variables are 
used. Rescaling the strain rate leads to the Weissenberg 
number Pe= 7r. Its values indeed indicate that both 
systems are driven with comparable rescaled rates, with 
the colloid experiment covering a broader but overlap- 
ping range in Pe. The similarity in the stress vs strain 
curves shown in Fig. [5] is recovered in the similarity of the 
parameters obtained when the schematic model is fitted 
to the data (Tab.[I|. Over the broader range in effective 
shear rates in the colloid experiment, slight drifts in the 
model parameters are visible which are not required for 
the metallic data plj . 

To reveal the microscopic origins of the stress-strain 
curves we investigate the local structural distortions. The 
universality in the metal and colloid rheology allows us 
to generalize findings in colloidal dispersions. Using con- 
focal microscopy we can observe individual particles and 
hence their rearrangements under flow, and the affine 
and non-affine motion can be separated. Figure [31 shows 
pair-correlation functions [29| determined at four strains 
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FIG. 2. Comparison of stress vs strain measurements for vary- 
ing strain rates in a metallic melt (a = axx, blue triangles) 
and a colloidal dispersion (a = a^y, red circles). The metal- 
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lie glass data is taken from Fig. [T] b) at T 
and the colloidal data of a fluid state is taken from Ref. [2l|]. 
Strain rates, converted to and listed as Weissenberg numbers, 
increase along the arrow. Solid lines are MCT fits with the 
schematic model for the parameters in Tab. [J 



after switch-on of shear in PMMA-coUoids (average ra- 
dius i?=780nm, polydispersity « 6%) slightly below their 
glass transition packing fraction (0 = 0.565 relative to 
4>c = 0.58 [lo| : same samples and confocal microscopy as 
in [3^). The shear distortion g(r) — g°^(r), i.e. the differ- 
ence between the sheared and equilibrium pair distribu- 
tion functions, exhibits an anisotropic quadrupolar and 
an isotropic angular contribution, which are obtained by 
angular averaging with the appropriate spherical func- 
tions [3H . The quadrupolar term relates to the shear 
stress and becomes more pronounced with strain; this is 
consistent with a very weak shear stress overshoot mea- 
sured in comparable PMMA systems [s^ . It describes a 
shift of the neighbor peak to smaller (larger) r along the 
comprcssional (extensional) axis. The isotropic distor- 
tion in g(r) — g'^'^{r) is unexpected from continuum elas- 
ticity, yet, for distances connected to the local caging, it is 
comparable to the quadrupolar distortion in magnitude. 
It corresponds to a shift of the neighbor peak towards 
smaller distances r, which we interpret as a local com- 
pression of the system. The very small compressibility 
of the system suggests a strong increase of the diago- 
nal parts of the stress tensor which give the pressure; 
it has been linked to shear banding jssl . [13], although 
flow remains homogeneous in our study. This could be 
caused by the different transient evolution (shear stress 
overshoot vs monotonous pressure increase). 

To support our finding of strong local compressional 
effects under shear we turn to simulations. A two- 
dimensional binary mixture of hard disks in a glassy state 
slightly above its MCT glass transition packing fraction 
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FIG. 3. Distortions of the pair correlation function, g[v) — 
g'"^{r) as a function of distance r: Quadrupolar 322 in panel 
(a), isotropic go — 17°'' in panel (b), for strains as labeled after 
switch-on of shear in PMMA colloids with average radius R; 
the shear rate corresponds to Pef^^ 79. 



pressure (not shown) , because of the very low compress- 
ibility of the system at the high packing fraction. Micro- 
scopic MCT calculations in two dimensions which extend 
the work in Ref. [s^ are included as inset in Fig. |4] and 
exhibit the non-monotonic behavior in the quadrupolar 
contribution but neglect an isotropic contribution [27 1 . 
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(0 = 0.81 relative to c^c = 0.80 [s^) exhibits stress-strain 
curves which share all the discussed features of the ex- 
perimental data [21| . Because the neighbor peak is fixed 



at the diameter d for hard spheres and cannot move in- 
ward, the observed compression should correspond to an 
increase in the contact value. It also determines all ele- 
ments of the stress tensor in the Virial-law: 

nkeTl - ^BTn'tf jdSrvgir). (4) 



J dS rr g{r) 



Here, g{r) accounts for the probability of collisions across 
the total internal contact area of the particles (s^l- The 
relation (j?]) lends itself to an efficient evaluation-scheme 
in event driven simulations [37| . To resolve the local con- 
tact values g{r), the collisions can be recorded to yield: 
5(f) - l/{2TTnd^NkBT) X (1/r,) Ecoii. |Ap(f)|, with the 
sum of the magnitudes of the momentum transfers in 
binary particle encounters in direction f in the time- 
interval Tc [ill. Figure U] shows the angle dependent 
contact value g{r) which is consistent with our experi- 
mental findings (data not shown); the angle starts from 
the flow-axis. Curves for different strains are shown. 
In the linear regime, where a^y = Goo 7s holds, the lo- 
cal distortion exhibits the expected quadrupolar angu- 
lar dependence. The contact probability decreases in 
the extcnsional and increases in the compressional di- 
rection. This behavior initially grows with strain. Be- 
yond the linear regime, the isotropic distortion develops 
and matches the quadrupolar dependence in magnitude; 
the latter varies non-monotonously and follows the stress- 
overshoot [31I, [s^l- In the stationary state, the contact 
probabilities are increased above the quiescent value for 
all angles. Even in the extensional direction, particles 
are pushed together by the shear flow. This is accompa- 
nied by a strong (monotonic) rise in the flow-dependent 



FIG. 4. Contact value g{r) as a function of angle 9 to the 
shear flow direction from Brownian dynamics simulations of 
a binary disk mixture in d = 2 at PeS> 1. An arrow indicates 
increasing strains as labeled; the black dotted line shows the 
equilibrium value. The inset shows corresponding MCT re- 
sults for hard disks where an isotropic term is neglected [27l ]. 

In conclusion, we have shown that bulk metallic glasses 
close to their glass transition temperature and colloidal 
dispersions close to their glass transition packing frac- 
tion exhibit equivalent stress-strain relations when prop- 
erly rescaled. Stresses agree when scaled by an entropic 
term characteristic for repulsive interactions; it is propor- 
tional to the thermal energy per particle volume. Fur- 
thermore, strains agree when geometrical factors aris- 
ing from the rheometer are accounted for, and shear 
rates agree when properly rescaled by the quiescent in- 
ternal relaxation time. Initially, stress increases linearly 
with strain, but for large Weissenberg numbers, when 
the external rate exceeds the inverse internal relaxation 
time, an overshoot originates due to negative stress auto- 
correlations. A schematic model of MCT rationalizes the 
rescaled curves with comparable parameters for metallic 
and colloidal systems. The microscopic origin of this ef- 
fect is investigated by microscopy, simulations and micro- 
scopic MCT. The quadrupolar distortion expected from 
continuum elasticity theory, dominates only during the 
linear regime. In the nonlinear regime, a strong rise in 
an isotropic distortion arises from the constraints of lo- 
cal packing. It corresponds to a local compression under 
flow and causes a monotonous increase in the pressure. 
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